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Abstract — In this work we compare the predictions of two 
representative hadronic interaction models, EPOS 1.99, and 
QGSJET 11-03 with several extensive air showers (EAS) pa- 
■ rameters for proton and iron primaries in the energy range 
10i7 . iQ^eV using CORSIKA-6990. The EAS parame- 
ters depth of shower maximum, shower size, size of muon 
shower, muon number distribution, electron number distri- 
bution, size of hadron shower, hadron energy sum, electron 
muon correlations, and, hadron energy spectra are studied 
in this paper. 

Index Terms — hadronic interaction models, EAS, Xmax 

I. Introduction 

WHEN high-energy cosmic rays enter into the Earth's 
atmosphere they initiate cascades of secondary par- 
ticles producing extensive air showers(EAS). Information 
regarding the shower generating primary particle have to 
be derived from the registered information of secondary 
particles at observation level. The interpretation of prop- 
erties of primary radiation derived from air shower mea- 
surements depends on the understanding of the complex 
processes of high-energy interactions during the develop- 
ment of air showers. From the number and distribution of 
various ground particles of the EAS, the reconstruction of 
energy and the mass of the primary particle can be done. 
' But to relate the observables to primary energy and mass, 
more reliable algorithms and detailed air shower simula- 
' tions are needed. By comparing the predictions from sim- 
ulation with measurements one can draw conclusions on 
I the primary mass composition of the arriving particles. 

Again predictions from simulations suffer from system- 
[ atic uncertainties mainly due to statistical fluctuations in- 
. volved in large-scale experiments and due to modeling of 
' HE interactions. While the electro-weak interaction pro- 
cesses are reasonably well understood; above the attained 
energy by the man-made accelerator, modeling of hadronic 
multi-particle production is subject to large theoretical un- 
certainties. Estimation of these uncertainties are further 
difficult task. 

Different hadronic interaction models predict different 
lateral shapes and different number of particles at obser- 
vation level. Hence, it is possible to test and compare the 
available interaction models by studying these EAS param- 
eters. Moreover, some of these parameters like the muon 
content, depth of shower maximum Xmax , the expected 
lateral shape etc. also depend on the mass of primary cos- 
mic rays. Heavier primaries lead, on average, to a flatter 
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distribution and lower value of Xmax and more muons. 

II. EPOS 1.99 AND QGSJET 11-03 

CORSIKA [1], [2 is a detailed Monte Carlo program 
to simulate the 4-Dimensional evolution of EAS in the 
atmosphere initiated by hadron, photon or any other 
particle. In CORSIKA-6990 ,3., there are seven HE 
models, namely DPMJET 2.55 ' @], EPOS 1.99 [5], [6], 
neXus 3.97 12], QGSJET-OIC [1], QGSJETH.3 [9], [10], 
SIBYLL 2.1 [n], [H], and, VENUS 4.12 [13], from which 
we can choose one model at a time for EAS simulation. 

To explain hadron-hadron/nucleus- nucleus collisions 
above the energy attainable by the man made accelerators, 
reliable and consistent hadronic interaction models are to 
be adopted. DPMJET, EPOS, QGSJET, and, SIBYLL 
are based on the Gribov's Reggeon approach [T^, [l^, of 
Pomeron exchange in multiple scattering. The Pomeron 
corresponds to microscopic partons (quarks & gluons) cas- 
cades, can be classified into "soft" , "semi-hard" , and 
"hard" Pomerons. Soft non-perturbative interactions in- 
volving large impact parameters and slow energy rise are 
described as "soft" Pomeron exchanges, are dominant at 
relatively low energy domain giving important contribu- 
tions to total, inelastic, and diffractive cross sections. How- 
ever, in the high energy regime "semi-hard" Pomeron ex- 
change, and, "hard" Pomeron exchange describing hard in- 
teractions with high energy rise, at large impact parameter, 
and at small impact parameters respectively are dominant. 

A. QGSJET 11-03 

QGSJET [8] is an addition of Gribovs Reggeon ap- 
proach [14], [15] of hadronic and nuclear collisions to 
the Quark-Gluon String model of high energy interac- 
tions. It has been generalized to treat nucleus-nucleus 
interactions and "semi-hard" processes using the so-called 
"semi-hard" Pomeron approach. , [10] , [16] . The "semi- 
hard" processes can be described by enhanced Pomeron 
diagrams and proved to be of extreme importance for a 
correct treatment of very high energy hadronic interac- 
tions [1^. QGSJET scheme is based on the assump- 
tion that individual Pomeron exchanges occur independent 
of each other, which is not true at high energy regime. 
At high energies parton cascades strongly overlap and in- 
teract with each other. These effects can be described 
as Pomeron-Pomeron interactions. This non-linear inter- 
action effects are incorporated into the newly developed 
hadronic interaction model QGSJET-II-3; which is based 
on the assumption that corresponding effects are domi- 
nated by "soft" partonic processes. Re summation of es- 
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sential enhanced contributions corresponding to particular 
final states of the interaction from uncut diagrams (repre- 
senting the elastic scattering), and from various unitarity 
cuts of enhanced Pomeron diagrams for all orders, yields 
the final state. QGSJET 11-03 is based on the obtained 
solutions, explicitly treating the corresponding effects in 
individual hadronic (nuclear) collisions [5], [TU] . 

B. EPOS 1.99 

EPOS [S]is a newly developed model emerging from 
VENUS [13] and neXus models [7]. EPOS is a par- 
ton model, with many binary parton-parton interactions, 
each one creating a parton ladder. EPOS is a consistent 
quantum mechanical multiple scattering approach based 
on partons and strings, where cross sections and the parti- 
cle production are calculated consistently, taking into ac- 
count energy conservation in both cases (unlike other mod- 
els where energy conservation is not considered for cross 
section calculations) [TS]. In EPOS 1.99 reduction of the 
proton-nucleus cross section is done for better correlation 
between the number of muons and the number of electrons 
at ground based air shower measurement. [19], [20] 

C. FLUKA 2011 

In CORSIKA apart from seven high energy interac- 
tion models there are three low energy models, namely 
GHEISHA 21 ,UrQMD [22j, and FLUKA [23 • The Fluka 
hadron-nucleon interaction models are based on resonance 
production and decay below a few GeV , and on the Dual 
Parton model above. Two models are used also in hadron- 
nucleus interactions [211, [Ml- 
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Average electron shower size Vs log of primary- energy 



D. Moun Number 

Muons are produced mainly by decay of charged pions 
and kaons in a wide energy range. Usually they are not 
produced directly on the shower axis. Multiple Coulomb 
scattering in the atmosphere and in the detector shielding 
may change the muon direction. It is known that the re- 
construction of the longitudinal development of the muon 
component provides a powerful tool for primary mass mea- 
surement, giving an information similar to that obtained 
with the fiuorescence technique, but in the energy range 
above that accessible by the detection of fluorescence light. 
In some experiment like KASCADE truncated muon size is 
calculated by integrating muons between 40m and 200m. 
Instead of total number of muons, truncated muon size is 
considered for E AS study [26] . The range of muon trunca- 
tion is from 140m-360m in KASCADE Grande. 
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Average truncated muon number Vs log of primary energy 
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E. Depth of Shower max 

The depth of shower maximum contains information 
about the mass of the primary CR initiating the shower 
as weh as about the properties of hadronic interactions. 
The average value X^ax depends on the primary energy E 
and on the number of nucleons A: 

X„,ax = a {InE- In A) + 13 (1) 
with a and (3 depending on the the details of hadronic 
interactions. Their values are very sensitive to changes in 
cross-section, multiplicity and elasticity |571. Eq.l can be 
derived from the simple generalized Hcitler model to show- 
ers induced by hadronic primaries, but it is in good agree- 
ment to the description of the Xmax evolution predicted 
by hadronic models currently in use. In some specific en- 
ergy range, a and /? may be considered independent of 
E. Another sensitive parameter is cr^^^^ , expressing quan- 
titatively the shower to shower fluctuations of X^ax ■ It 
depends mainly on the cross section and less strongly on 
the elasticity. This makes fluctuations in Xmax , a good 
parameter to study hadronic cross sections at ultra-high 
energies. p7] 

III. Simulation 

The shower simulations arc performed using CORSIKA- 
6990. Hadronic interactions at low energies(E < 80 GeV) 
are modelled using the FLUKA-2011 code. High Energy 
interactions are treated with EPOS 1.99 and QGSJET H- 
3. Vertical showers initiated by primary protons and iron 
nuclei are simulated. Observation level is taken at the 
sea level, US standard atmosphere taking default mag- 
netic field in CORSIKA. 1000 showers are simulated each 
for three primary energies 10^^ , 10^^ felO^^e^, two pri- 
mary mass(p & Fe) and for two HE models EPOS 1.99 and 
QGSJET H-3. All together 12,000 showers are generated. 
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IV. Results 

The Depth of maximum Xmax , and the RMS value 
'^Xmax f'^'' the two models are plotted as function of en- 
ergy in Figures [H It is seen that heavier primary produces 
shower maximum at lower depth with lower fluctuation 
compared with lighter primary as expected. The aver- 
age total number of electrons (A^'e), muons(A^p), truncated 
muons(A^^''), hadrons {Eh > lOOGey), and sum of energies 
of all the hadrons(£^/i > lOOGeV) registered at the ground 
level for the interaction models EPOS 1.99 and QGSJET 
11-03 are plotted as a function of energy in Figures [21 El HI 
[3 and m respectively. 

It is seen that, both models yield a linear dependence of 
these components with energy in Log-Log scale. 

A. Dependence of Xmax with primary Energy 

The value of Xmax and also the slope of the curve for 
proton primary for EPOS 1.99 is slightly higher than that 
of QGSJET 11-03. Again, for iron primary also the slope 
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Distribution of the number, of truncated muons for proton 
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TABLE I 
Slope of Xmax^^o Intercept 



Model 


Slope of Xmax 


Intercept 


& Primary 


a 


error 


P 


error 


EPOS-p 


58.3690 


±1.64891 


-303.488 


±29.7108 


QGSII-p 


47.5875 


±0.24876 


-125.219 


±4.48227 


EPOS-Fe 


62.9539 


±1.55497 


-491.491 


±28.0169 


QGSII - Fe 


55.2977 


±0.23593 


-354.504 


±4.25079 



of the curve has higher value for EPOS 1.99, while the 
values of Xmax are nearly equal for both the models. The 
value of Ux^^^ for proton primary for EPOS 1.99 is slightly 
higher(average about 3% high) than that of QGSJET II- 
03. But for iron primary there is no such trend observed. 
The fitted values of the slope a and intercept (3 are given 
in the Table HI 



B. Primary Energy Correlations 

From Figures [5] through ini it is seen that, there is no sig- 
nificant differences in the electron shower size and sum of 
hadron energy for the two models considered in the chosen 

TABLE II 

Average Truncated Muon Number and RMS value 



Energy/Model 


Proton Primary 


Iron Primary 


(eV) 




< «s(Jvr) > 


RMS 


< hiK') > 


RMS 


IQl' - 


QGSII 


5.33582 


0.09227 


5.45870 


0.04113 


IQi^- 


EPOS 


5.42664 


0.10664 


5.52427 


0.04210 


1018- 


QGSII 


6.26616 


0.09035 


6.38253 


0.04637 


1018- 


EPOS 


6.37158 


0.09563 


6.47433 


0.04434 


10"- 


QGSII 


7.21209 


0.08441 


7.31510 


0.05116 


10"- 


EPOS 


7.31632 


0.09665 


7.41770 


0.04526 



TABLE III 

Average Electron Number and RMS value 



Energy/Model 


Proton Primary 


Iron Primary 


(eV) 




<ig(Ne)> 


RMS 


< ig{Ne} > 


RMS 


lOl'- 


QGSII 


7.28989 


0.14648 


7.05567 


0.08888 


1017- 


EPOS 


7.27924 


0.15111 


6.99535 


0.09171 


1018- 
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8.40937 
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Fig. 8 

Distribution of the number of truncated muons for iron 
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Distribution of the number of electrons for proton primary 



primary energies for proton and iron primaries. However 
for muons (Fig. [31 S]) differences between the predictions 
of the two models are significant. The number of muons is 
larger for EPOS1.99 model as compared to QGSJET 11-03 
for both the primaries. From the figure [3l it is seen that 
there is very little overlap in the region bounded by p & 
Fe primaries for two models considered. However, for trun- 
cated muons (in the range 40m-200TO) there is no overlap 
beyond 10^*^ (Figure S]). For Eh > lOOGeV hadrons, to- 
tal number for both proton & iron primaries are found to 
be significantly more for EPOS than that for QGSJET II-3 
(Figure EJ. However, considering the sum total of energies 
of all the hadrons {Eh > lOOGeV), there is no significant 
difference between the two models (Figure [B]). 



C. Distribution of truncated muon numbers and electron 
numbers 

Distributions of truncated muon numbers lg{N*^) for 
proton and iron primaries are plotted in Figure [7] and [8] It 
is seen that EPOS1.99 yields higher value of lg{N^J^) than 
QGSJET 11-03. The mean & s.d. of these distributions 
are tabulated in Table [Til and significance test done (|C|) . 

Distributions of electron numbers lg{Ne) for proton and 
iron primaries are plotted in Figurel9landfT0l It is seen that 
EPOS1.99 produces slightly less numbers of electrons for 
energy lO^el^for iron primary, but at lO^^eV^, EPOS1.99 
produces slightly more electrons as tabulated in Table 1111) 
and significance test done([Cl). 
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Distribution of the number of electrons for iron primary 
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Number of electrons Vs number of truncated muons 



D. Electron to muon correlation 

The average number of electrons as a function of the 
number of truncated muons are plotted in figure [11] for 
the two models. Slopes of all the curves are almost 
equal. To emphasize the differences between the model 
predictions, relative deviation in the model prediction 
EPOS1.99 with respect to the QGSJET II-3 prediction, 
{N^^^^ - N^^^)/N^^^ is plotted against the mean of 
j^QGS ^^^^^ j^EPOS jj.^ Figure [H It is seen that for pro- 
ton primary, EPOS1.99 predicts slightly less (about 2%) 
electrons for lO^^'et^ proton-induced showers, slightly more 
(about 3%) electrons for lO^^eV and about 8% more elec- 
trons for lO^^eV primary energy. But for iron primary 
at IQi^eF, lO^^eV, and, lO^^eV] EPOS 1.99 yields about 
12% lesser, about 6% lesser and nearly equal ( difference 
is less than 1%) numbers of electrons as compared with 
QGSJET II-3 predictions. 

E. Variation in hadron component 

Total average number hadrons with energy Eh > 
lOOGeV, for both primary masses & at all the three en- 
ergies show significant difference between the two models 
(|D|) . Energy spectrum of registered hadrons for 10^"^ eV, 
lO^^eV, and lO^^eV primary for both the models consid- 
ered are displayed in the Figures [T31 HH and [13] respec- 
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Relative deviation of the number of electrons Vs average 
number of truncated muons 
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HADRON ENERGY DISTRIBUTION (10-*- '^eV) 



tively. There are apparently no distinguishable differences. 
All the four data plots are overlapped to each-other with 
their error-bars. 



V. Summary 

Although Xmax , shower size, hadron energy sum, 
hadron energy distribution shows no significant difference 
between the two HE models EPOS 1.99 (FLUKA) and 
QGSJET II-3 (FLUKA); muon number, hence electron 
to muon correlation shows incompatibility between them. 
EPOS 1.99 (FLUKA) predicts more muons than QGSJET 
n-3 (FLUKA). Also hadTon{Eh > lOOGeV) shower size 
prediction by EPOS 1.99 (FLUKA) yields higher values 
irrespective of primaries considered as compared to that 
by QGSJET n-3 (FLUKA). Modifications may be needed 
for one, or, both of the models considered herein by com- 
paring the model predictions with experimental data. 
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A. Statistical Test 

A. Z-test results for X. max distribution fTable [TV\) 

The null hypothesis is 
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(QGSJET 11-03 data and EPOS1.99 data). 
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Xmax Dist 



Inference 



Not 

Moderately 

Not 
Highly 
Moderately 



gnificant 
gnificant 
gnificant 
gnificant 
gnificant 



Significant 



B. Z-test results for a and j3 (Table VV\) 

The null hypothesis is 
iJo = There is no difference between the parameters for 
QGSJET 11-03 and EPOS1.99 . 
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TABLE V 

Z STATISTICS FOR a AND /3 



Parameter 


a and /3 


k^Primary 


Z 


Inference 


a — Proton 


6.466 


Significant 


a — Iron 


4.981 


Significant 


P — Proton 


5.933 


Significant 


P — Iron 


4.834 


Significant 



It is seen that the parameters a and j3 describing 
Xmax (Equation 1), show significant differences for EPOS 
and QGSJET II, irrespective of the primaries. 

C. Z-test results for N^^ and distributions (Table VvJ^) 

The null hypothesis is 
-ffo — There is no difference between the two samples 
(QGSJET n-03 data and EPOS1.99 data). 



TABLE VI 

Z STATISTICS FOR, A^^'" AND A^e DISt'^S 



Energy/Primary 
{eV) 


N" 




Nc 


Dist"- 




Inference 


7 


Inference 




— Proton 


20.36 


Highly 


1.60 


Not 








Significant 




Significant 




— Iron 


35.22 


Highly 


14.93 


Highly 








Significant 




Significant 


1018 


~ Proton 


25.34 


Highly 


2.86 


Moderately 








Significant 




Significant 


1018 


— Iron 


45.25 


Highly 


7.56 


Significant 








Significant 






1019 


— Proton 


25.68 


Highly 


7.93 


Significant 








Significant 






10" 


— Iron 


47.50 


Highly 


0.56 


Not 








Significant 




Significant 



Average truncated muon number shows highly signif- 
icant difference between the two models. For the aver- 
age electron number, for lower energy(lO^^ey) high sig- 
nificance for iron primary, but no significant difference for 
proton primary is seen, whereas at higher energy(lO^^ey), 
no significance for iron, but significant differences for pro- 
ton primary is observed ( Table IVI[) . 



D. Z-test results for total hadron number (Eh > lOOGeV) 
(Table\Vl^ 

The null hypothesis is 
Hq — There is no difference between the two samples 
(QGSJET 11-03 data and EPOS1.99 data). 

TABLE VII 



Z statistics for total hadron number (Eft > lOOGeV) 



Energy/Primary 


Nh Disf^ 


ieV) 


Z 


Inference 


1017 _ Proton 


4.07 




Si? 


>nificant 


10^^ — Iron 


21.09 


Highly 


Si? 


jnificant 


IQis _ Proton 


8.91 


Highly 


Si? 


>nificant 


IQis _ Iron 


6.52 


Highly 


Si? 


>nificant 


1019 _ Proton 


13.25 


Highly 


Si? 


jnificant 


10^^ — Iron 


15.59 


Highly 


Si? 


>nificant 



